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Abstract. After discussing some hints for possible masses of quasiparticles in 
quark matter on the basis of lattice equation of state, we present pion and kaon 
transverse spectra obtained by recombining quarks with distributed mass and 
thermal cut power-law momenta as well as fragmenting by NLO pQCD with 
intrinsic kj* and nuclear broadening. 



1. Hadronization by quark coalescence 

Quark coalescence or recombination proves to be a remarkably simple and powerful 
model of hadronization of quark matter in central, relativistic heavy ion collisions. 
Transverse spectra and the constituent quark scaling of ellipticity in the collective 
flow are predicted well by this concept, although theoretical arguments have been 
raised against it time to time. In this contribution we review first a few arguments in 
favor of the quark recombination model ALCOR and then present application results 
to pion and kaon transverse spectra at RHIC. We propose to consider the sliding slope 
of transverse spectra, with transverse flow correction wherever it seems to be necessary, 
as an excellent indicator of deviations from exponential and as a smooth fitting edge to 
high-p*r spectra calculated in NLO pQCD by the use of fragmentation functions. This 
sliding slope is linear in the energy in a co- flowing cell for a cut power-law distribution, 
so it reflects the power-law tail nicely. The inverse powers, the l/(q — 1) parameters, 
on the other hand should scale according to quark coalescence ratios for mesons and 
baryons when produced by recombination. 

A common worry with respect to recombination is represented by the question 
of entropy reduction due to diminishing the color degree of freedom. Albeit the 
dynamical details of any non-perturbative hadronization process are unknown to date, 
there can be constraints from lattice QCD equation of state (eos) on the possible 
reduction. In ref.pQ we have shown that lattice QCD allows for a significant reduction 
in the effective number of ideal gas particles by a constant total entropy (adiabatic 
process). This way the recombination process does not necessarily reduce entropy, if 
at the same time an adiabatic expansion of the volume takes place. 

Another point is raised by asking why to bother statistical and thermal models 
of quark matter when perturbative QCD (pQCD) is able to predict both the equation 
of state and the pr spectra. The grounds here are also theoretical: on the one 
hand non-perturbative methods in field theory are available for thermal systems to 
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Figure 1. The distribution of invariant mass squared for two Boltzmann 
distributed, massless partons. As examples an IR safe quantity, M 4 , and an 
unsafe one, l/ln(M 2 ) are also indicated by thinner lines. 

date, therefore it is important to investigate whether the phenomena we observe 
experimentally may be thermalized or not, on the other hand pQCD might not be 
applicable for infrared unsafe physical quantities in a thermal state at all. 

To illustrate this, figure Q] plots the probability distribution function of the 
invariant mass square, M 2 , combined from two massless, Boltzmann distributed 
partons: 

P(M 2 ) = N J dVd 3 ^ h h 5 (M 2 - 2E 1 E 2 + 2pi • p 2 ) , (1) 

with pi = (Ei,pi) four-momenta of the partons, fa = e~ Ei l T Boltzmann distributions 
and a temperature dependent normalization factor M . The probability is normalized 
as J P(M 2 )dM 2 = 1 and due to kinematic reasons < M 2 < AE\E 2 for physical 
particles. The result of the integral ([J) is given by 

pc^ w 

This probability is finite at low M, a typical infrared unsafe quantity is l/ln(M 2 /A 2 ). 

For the hadronization from a quark matter we shall consider cut power-law 
(Tsallis) distributed quark momenta at a given temperature, T, and parameter q: 

W)4(l + (,-D^)^. 0) 

The energy - momentum dispersion relation is presumably not near to the free one 
at the hadronization process. In the distributed mass model we consider a massive 
dispersion relation, E(p) = y m 2 + P 2 and a mass distribution, p(m). The smeared 
quark masses give some flexibility to recombine any hadron mass we target at, in 
particular the low pion mass. We use two approaches: i) a model mass distribution 
suppressing both low and high quark masses [3], 
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Figure 2. Interaction measure scaled by T 2 for pure SU(3) lattice gauge theory 
[5] and for 2+1 flavor QCD with physical pion mass obtained recently [§]. The 
constant high-T asymptotic of this quantity may correspond to a constant mass 
[4] or to screened strings contributing to the free energy density as n 2 / 3 [7|. 



and ii) a distribution with a mass gap deviced to mimic lattice eos data [5] (see later) , 

, , 4 TOn / m 2 . , , 

p(m) = ^V 1 -^- 5 

Before presenting our results for the meson spectra, we discuss briefly what lattice eos 
may hint at the existence and shape of p{m) distributions. 



2. Quasiparticles and QCD eos 

We would like to point out that T > T c quark matter does not necessarily behave like 
approaching an ideal quark-gluon plasma (QGP) at very high temperature. The best 
test considered for this is the interaction measure, e — 3p, which vanishes for an ideal 
massless gas and in the T — > oo limit. Lattice gauge theory and QCD data all agree 
with this expectation. 

However, a different question is the high-temperature scaling of this quantity. 

Figure [2] plots the interaction measure scaled by the square of the temperature 
(not, as usual, by T 4 ) for pure SU(3) lattice eos [5] and for a more recent 2+1 flavor 
QCD eos with physical pion mass [6]. Both seem to agree with a constant for 




soon after the color deconfinement up to the highest temperature investigated on the 
lattice. 

Not any model of quark matter satisfies this constraint. In the quasiparticle 
picture for an ideal gas of mass m particles the T — > oo limit is just m 2 /2, in case 
of a mass distribution it is given by (m 2 )/2 — d/dT((m 2 ))/(4T 5 ) to leading order. 
While this may well be a constant for some temperature dependent mass distribution 
p(m; T), the pQCD leading term, T 3 d/dT(l - a 2 g 2 ) = (a^/ao) g 2 T 2 for the usual one 
loop expression g 2 (T) = oto/ ln(T 2 /A 2 ), is not a constant 

| although also proportional to the thermal mass squared. 
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Figure 3. Energy distribution in a blast wave's rest frame for pion pt spectra 
(upper) and the sliding slope obtained from it (lower). 



In principle in advanced quasiparticle models the mass smearing may also be 
temperature dependent. In this case a mean field term occurs in the pressure, which 
nevertheless cancels in the combination e + p. Integrals of the PDF p{m) leading to 
a given e + p can be compared to lattice eos, too. In this case the Markov inequality 
gives estimates for the lowest mass [J. 

3. Hadron pt spectra from quark matter 

Finally we present transverse momentum spectra for pions and kaons as obtained from 
quark recombination and pQCD fragmentation. Figure [3] and figure [?] show for pions 
and kaons the respective curves. The upper patterns are invariant yields, f(E), while 
the lower pictures are the sliding slopes 

Aslope = - (^W(£)) • ( 7 ) 

Here E = r y(mT — vpx) is the average transverse flow corrected energy at zero rapidity. 
We took into account a transverse flow of v w 0.5 for the kaons. For a Tsallis 
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distribution T s \ ope (E) = T* + (q — l)E, the steepness of the linear rise is not affected 
by the blue shift factor, but the intersect is: T* = Ty/(l + v)/(l - v) for E > m. 

We used a next-to-leading-order pQCD-improved parton model for calculating the 
high-px pion and kaon spectra in AuAu collisions. Theoretical calculations presented 
here, are based on 'standard' parton distribution functions [10] (PDF), but in a 
generalized way: a 2-dimensional Gaussian transverse momentum distribution were 
included, represented by the average intrinsic transverse momenta, (k^) (see detail 
in Refs. [3 [9j). The calculation were carried out both without (standard pQCD 
calculation) and with intrinsic- fey, as motivated by experimental data [111112] . 

The calculated spectra strongly depend on initial and final state nuclear effects, 
which have to be taken into account in proton-nucleus and nucleus-nucleus collisions. 
Here we used the nuclear shadowing and multiple scattering, and the jet-quenching 
with the corresponding opacity parameter value, L/X » 3.5 — 4.0. 

We applied two parameterizations for fragmentation: KKP [T3] and the latest 
version of AKK [T4] , 

In conclusion the curves from thermal cut power-law distributed quark matter 
recombination meet smoothly those NLO pQCD results, which take into account an 
intrinsic transverse momentum in the parton distributions. The (q — 1) steepness of 
the energy dependence of the inverse logarithmic slope, T s iopc is expected to follow a 
coalescence rule for mesonic and baryonic recombination at high px only. 
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Figure 4. Conjectured co-moving energy distribution from STAR[T5] kaon pt 
data (upper) and the sliding slope calculated from it (lower). Theoretical points 
from quark matter recombination and lines from NLO pQCD fragmentation 
calculation are plotted. The lower line is without intrinsic fcy assumption, for 
the upper line we assumed (fcJ,) = 4 GeV 2 /c 2 . 



